Electrothermodynamic power converter with converging flows

ABSTRACT

An electrothermodynamic generator and a compressor is described for use in a Marks/Ericsson Cycle. It uses a charged liquid tin droplet/nitrogen aerosol for heat/electric conversion at a maximum temperature of about 1800 K and a charged water droplet/nitrogen aerosol electrothermodynamic compressor at 300 K in a Marks/Ericsson Cycle which has a theoretical efficiency of 83%, and a practical efficiency of 60% to 70%. An ETD compressor with a water charged droplet is described. The other cycle components, comprising a heat source, liquid tin loop, heat exchangers, water loop and liquid pump are not here considered.

1. FIELD OF THE INVENTION

This invention relates to an electrothermodynamic generator with converging gas and charged aerosol flows. Compared to prior art devices, the present invention enables the direct conversion of the thermal power of a gas aerosol to electric power, at several orders of greater magnitude of electric power and current, at about the same voltage, and in the same volume.

2. BACKGROUND OF THE INVENTION

A comprehensive background and bibliography in the field of this invention has previously been given.¹ This invention results from further studies in connection with Methods III.²,2.1 and IV.¹

3. IN THE FIGURES

FIG. 1 shows a cross section of an electrothermodynamic generator of this invention and a diagram of the electrical connections, showing the outer curve for the surface of the converging duct, and the inner curve for the surface of the charged aerosol within the flowing gas in the duct.

FIG. 2 shows a graph of velocity ratio (U_(z) /U_(b)) vs. distance ratio (z/b)

FIG. 3 shows a cross sectional view of an electrothermodynamic generator of this invention comprising a 9 element array.

4. DESCRIPTION OF THE INVENTION

In FIG. 1, there is shown a straight duct 1, which contains a flowing gas 2 which passes into a converging duct 3, and thence through an orifice 4. A small diameter tube 5 within the duct 1 contains a flowing conductive liquid 6 such as molten tin, Wood's metal, mercury, water, an alcohol, etc. which issues as a liquid jet 7. An exciter electrode 8 is located at the entrance to the converging duct 3. A voltage source 9 is connected across the exciter electrode, and the tube 5, applying an electric field between the jet 7 and the exciter electrode 8. As a consequence the liquid jet breaks into a large number of small diameter charged liquid droplets, thus forming a charged aerosol at the entrance to the converging duct 3, which has a diameter d_(a) less than the diameter D_(a), at z=a. The converging duct 3 converges along a first or outer surface along an outer curve 12. The charged aerosol expands to a smaller diameter d_(a) at z=a, and converges along a second or inner surface, which nests within the said first surface. The converging gas stream line has a velocity vector OA, which is resolved into an axial vector velocity component OB, and a normal vector velocity component AB directed towards the z axis. The normal velocity vector AB causes the charged aerosol to be centrally confined, and to be pinched around the axis within the duct; and to remain separated from the wall of the duct 3. This prevents the wall from being wetted by the liquid, and prevents "shorting" the charged aerosol potential difference via the walls of the duct. By this means, the gas and the charged aerosol are "uncoupled". The gas flow is consistent with the continuity laws, the charged aerosol is subjected to the velocity increase along the z axis provided by the flowing gas, but is contained within a different or inner surface envelope suitable for maximum power conversion. A further discussion of the invention will follow the Mathematical Physics Section.

                                      TABLE OF SYMBOLS                             __________________________________________________________________________                    Description              Units                                  __________________________________________________________________________     Symbol                                                                         A              Area                     m.sup.2                                C              Specific Heat            joules/kg                              D              Diameter of Duct         m                                      E              Electric Field Intensity v/m                                    I              Current                  amps                                   J              Exponent of Area for the Duct Surface                           K              A Constant                                                      L              Length Along Conversion Space = b - a                                                                   m                                      M              Mach No                                                         N              Number                                                          P              Power                    watts                                  R              Gas Constants - see constants, value of ;                                                               a ratio                                S              Entropy                                                         T              Temperature              °K.                             U              Velocity                 m/s                                    V              Potential, electric      volts                                  a              Distance from origin at z = a                                                                           m                                      b              Distance from origin at z = b                                                                           m                                      d              Diameter of Charged Aerosol                                                                             m                                      e              Electron Charge see constants, values of                                                                c                                      f              A factor                                                        j              Exponent of (z/b) for cross sectional                                          area of the charged aerosol surface                             k              A constant                                                      m              Mass                     Kg                                     n              Exponent of (z/b) for velocity                                  p              Pressure                 N/m.sup.2                              r              Radius                                                          z              A distance along the flow axis                                                                          m                                                     from the origin                                                 ν           Volume                   m.sup.3                                α        •m.sub.o /.m.sub.g ratio of liquid droplet mass flow                     rate                                                                           to the gas mass flow rate                                       δ        Density                  Kg/m.sup.3                             δ.sub.r  Relative density, a ratio δ/δ.sub.1 at                             273° K.                                                  γ        Ratio of specific heats c.sub.p /c.sub.v for gas                ρ          Electric charge density  c/m.sup.3                              η          Efficiency                                                      ζ         Electric field ratio (E/Ec) at incipient                                       breakdown                                                       Ω        Molecular weight                                                χ          Voltage ratio Ea/Ec                                             -χ         χ/ζ voltage ratio maximum                              .sup.--Γ χ/ζ.sup.2 power ratio maximum                          Δ        Increase or decrease of a value                                 SUBSCRIPTS                                                                     L              Length                                                          Q              Heat                                                            a              The distance z = a along conversion space z axis                               from origin                                                     b              The distance z = b along conversion space z axis                               from origin                                                     c              Referring to same variable computed for a constant                             area duct, constant velocity for electric field,                               voltage, current and power                                      d              Duct                                                            e              Electric or Electron                                            g              Gas                                                             k              Kinetic                                                         o              Droplet                                                         p              Pressure                                                        r              Relative value                                                  s              Sonic                                                           v              Volume                                                                         Air; at standard conditions, 273K and 1 atmos. =                               1.01 × 10.sup.5 N/m.sup.2                                 SUPERSCRIPTS                                                                   •        Time Rate                                                       --             Maximum value. When used thus -χ, .sup.--Γ these                     ratios                                                                         are the values at incipient electric breakdown                                 resulting from a η increase in the current by the                          ratio 1/ζ.                                                 CONSTANTS, VALUES OF:                                                          K.sub.e        1/2(ε.sub.o b.sub.1.sup.2)(Ω.sub.1 /δ.s                    ub.1)[1/273 × 8314)] = 4.00 × 10.sup.-4                            a constant of the electric power conversion                                    equation                                                        ε.sub.o                                                                               8.854 × 10.sup.-12 dielectric constant of                                                         f/mce                                  δ.sub.1  1.273 density of air at 273K and 1 atm                                                                  Kg/m.sup.3                             b.sub.1        3 × 10.sup.6 electric breakdown field at                                                          U/mK                                                  and 1 atmos., of air.                                           Ω.sub.1  28.97 average molecular weight of air                           R              8314/Ω gas constant                                                                               joules/Kg K°mol                 e              1.601 × 10.sup.-19 electric charge of the                                                         Clectron                               η.sub.k    1 - (U.sub.a /U.sub.b).sup.2 Kinetic power change in                           conversion                                                                     space as a fraction of the output kinetic                                      power                                                           __________________________________________________________________________

6. MATHEMATICAL PHYSICS 6.1 Differential Equation of the Conversion Space

The generalization of the earlier equations³ involved the use of a differential form of Gauss' equation in which the charge density p varies from z=a to z=b:

    d/dz(AdV/dz)=-Aρ/ε.sub.o                       (1)

in which the conversion space is between z=a and z=b, and

    A.sub.z ≡A.sub.a (z/a).sup.j                         (2)

    U.sub.z ≡U.sub.a (z/a).sup.n                         (3)

    p.sub.z =I/A.sub.z U.sub.z =(I/A.sub.a U.sub.a)(z/a).sup.-(n+j) =ρ.sub.a (z/a).sup.-(n+j)                             (4)

The differential equation obtained from (1)-(4) inclusive and to be solved is:

    d/dz[(z/a).sup.j (dV/dz)]=-ρ.sub.a (z/a).sup.-n        (5)

6.2 The Electric Field and the Current Ratio

The equation (5) is integrated subject to E=dV/dz=0 at z=b The equation for E is compared for various values of j and n to the Electric Field intensity E_(a) at z=a, the entrance to the conversion space, using the equation for E_(c) derived in the prior art⁴, for a constant velocity gas in a constant area duct.

    E.sub.c =-ρL/ε.sub.o =IL/ε.sub.o UA    (6)

    E.sub.b =b.sub.1 K.sub.a δ.sub.r =E.sub.c =IL/ε.sub.o UA (6.1)

An electric field intensity ratio ζ_(a) at z=a is derived: ##EQU1## For U increasing, n>1 and for (b/a)>>1, the numerator and equation (7) may be approximated: ##EQU2## and the field ratio ζ_(a) is small and large currents can be sustained without electric breakdown.

For the charged aerosol velocity decreasing n is negative η<0 and the approximation is ##EQU3## and here the field ratio ζ_(a) is large.

Hence it is preferred to use the conditions for (8). From (6.1) the maximum current I is limited by the breakdown field E_(b) =E_(c). According to Equation (7) the electric field intensity at z=a is decreased by an increased velocity gradient, that is for a positive value of η. Therefore according to (6.1) the current may be increased inversely as the ratio 1/ζ. Hence the maximum current ratio is 1/ζ, which, from (8) is

    1/ζ=(η-1)(b/a)                                    (10)

The value of n should be as large as possible but it is limited to attainable values of the ratio of the velocity at the exit plane to the velocity at the entrance plane of the conversion space; that is (U_(b) /U_(a))≈1000; for a useful ratio of the distance b to the distance a; that is (b/a)≈100. The useful maximum value of n is then from (3):

    η=ln (U.sub.b /U.sub.a)/ln (b/a)=1000/ln 100=1.5       (11)

From (3) and (10) the current ratio is:

    1/ζ=(η-1)(U.sub.b /U.sub.a).sup.1/η           (12)

For (U_(b) /U_(a))=1000 the maximum value of 1/ζ at η=1.2 is: 63.10. Putting n=1.2 and (U_(b) /U_(a))=1000 in (11), and solving for

    (b/a)=1η.sup.-1 (1/1.2) ln 1000=316                    (13)

The electric field intensity ratio ζ varies with (z/b) according to an equation derived from (5) by integration subject to the condition E=0 at (z/b)=1 ##EQU4##

This equation has a maximum value ζ for j=-η, but not for η+j≅1. If there is a maximum ratio ζ then its value and position is respectively given by (14) and (15): ##EQU5## and from (9) and (15) for (b/a)>>1

    (ζ/ζ.sub.a)=(-j).sup.-j (η-1)                (16.1)

A maximum exists only when ζ>ζ_(a), for which

    (-j).sup.-j (η-1)>1                                    (16.2)

For n=1.5 and (j+n)=0.84 in which case there is no maximum ζ.

However maxima do exist for cases where the continuity condition implies γ≈-η in which case

    (η-1)(-γ).sup.-γ =(η-1)η.sup.η >1 (16.3)

for positive n.

6.3 The Electric Voltage and Voltage Ratio

The voltage V from z=a to z=b is obtained by a second integration of equation (5), starting with equation (14), which is the result of the first integration. The second integration is subject to the condition that V=0 at z=a.

As a basis of comparison the voltage V_(c) across the distance L=b-a, for a constant area duct, with a charged aerosol at constant velocity, previously derived⁴, is used:

    V.sub.c =ΥL.sup.2 /2ε.sub.o =IL.sup.2 /2ε.sub.o Au (17)

A voltage ratio X is defined, and derived from the second integration: ##EQU6##

Equation (18) simplifies for (b/a)>>1 and η>1:

    χ=[1/(1-j)(2-η-j)](a/b).sup.η+j                (19)

From this, it is shown that, for η+j≈1, and j negative for a converging charged aerosol with increasing velocity, the voltage ratio is small.

However, from (17) the current I may be increased until incipient electric breakdown by the current ratio 1/ζ from equation (10); then χ≡χ/ζ is the corresponding maximum voltage ratio. Hence from (10) and (19)

    χ≡χ/ζ=[2(η-1)/(1-γ)(2-η-J)](b/a).sup.1-.eta.-γ                                                 (20)

Examples: ##EQU7## In the particular case where n+j=1

    χ=2(η-1)/η                                     (23)

For χ=1, equation (20) requires these values: ##EQU8## For χ=2, which may be about the maximum useful value: ##EQU9##

As shown in (12) a maximum current ratio is obtained 1/ζ=63.2 for η=1.2. However, by (23) χ=1/3; and so the power ratio Γ=63.2α(1/3)=21. If the voltage ratio χ increased to χ=1 by increasing the length by √3 then the power ratio Γ=36.4, not as gret as in (24).

For all values of χ, χ=1, χ=2, etc. The maximum current ratio 1/ζ is the same, since from (10) it depends only on the value of n. Thus if in (24) the voltage V=50,000 volts, in (25) the voltage is V=100,000 volts; the current being the same for each; but twice the power for (25), compared to (24).

6.4 The Electric Power and the Power Ratio

At incipient electric breakdown, the maximum voltage V across the converter from z=a to z=b occurs when the current is increased by the ratio 1/ζ, compared to the maximum current for the constant area, constant velocity case (6.1). The maximum power ratio is:

    γ=IV/I.sub.c V.sub.c =χ/ζ.sup.2 =χ/ζ(26)

If χ<1, the length L=(b-a) in (17) may be increased by χ^(-1/2). This decreases the current by χ^(1/2) increases the voltage to χ=1, and increases the power by χ^(-1/2) ##EQU10##

According to (11), the maximum practical values of (b/a)=100 and n=1.5; hence the maximum power ratio for χ=1 is from (23) and (12)

    Γ=1(η-1)(b/a)=1(1.5-1)100=50                     (28)

and for

    χ=2

    Γ=100                                                (29)

6.5 The Duct Surface

The duct surface contains the flowing gas which converges toward an orifice of diameter D_(b). Applying the continuity equation:

    mg=δ.sub.a U.sub.a A.sub.a =δ.sub.b U.sub.b A.sub.b =δ.sub.z U.sub.z A.sub.z                            (30)

Substituting (2) and (3) in (30) and solving for (δ_(a) /δ_(z))

    (δ.sub.a /δ.sub.z)=(U.sub.z /U.sub.a)(A.sub.z /A.sub.a)=(z/a).sup.η ·(z/a).sup.J =(z/a).sup.η+J (31)

when z=b

    (δ.sub.a /δ.sub.b)=(z/b).sup.η+J           (32)

The equation of the duct surface curve in the YZ plane is derived from (2). The j is capitalized J to distinguish it from the j used with the charged aerosol surface

    (D.sub.z /D.sub.a).sup.2 =(A.sub.z /A.sub.a)=(z/a).sup.J =(b/a).sup.J ·(z/b).sup.J                                     (33)

    D.sub.z =[D.sub.a (b/a).sup.J/2 ](z/b).sup.J/2             (34)

When (z/b)=1

    D.sub.b =D.sub.a (b/a).sup.J/2                             (35)

and since

r_(z) =1/2D_(z)

The equation of the duct curve is:

    r.sub.z =r.sub.b (z/b).sup.J/2                             (36)

The velocity U_(z) at (z/b) is from (3)

    U.sub.z =U.sub.a (z/a).sup.η =[U.sub.a (b/a).sup.η ](z/b).sup.η(37)

The exit velocity is U_(b) at (z/b)=1; hence

    U.sub.z =U.sub.b (z/b).sup.η                           (38)

The velocity space gradient at (z/b) is

    (dU.sub.z /dz)=η(U.sub.b /b.sup.η)z.sup.η-1 =η(U.sub.b /b)(z/b).sup.η-1                                      (39)

Example ##EQU11##

Hence the velocity gradient increases by a factor of 10 from z=a to z=b for this example.

6.6 The Charged Aerosol Surface

If the charged aerosol is uniformly distributed across the duct the continuity equation (30) must be used.

Example: Given δ_(a) /δ_(b) =1.5; b/a=100 Find (J+n) from (30)

    J+η=ln (δ.sub.a /δ.sub.b)/ln (b/a)=0.405/4.605=0.09 (40)

and if n=1.5 then J=-1.41. Putting these value in (20) χ=11.5 which is too large; hence the surface of the charged aerosol must be separated from the duct and j+n≈1 , so that χ=1 or 2.

When the surfaces of the duct and the charged aerosol are separated as shown in FIG. 1, the gas stream lines between the duct surface and the charged aerosol surface may be represented by a velocity vector, such as OA from the point 15. This vector has an axial component, and a normal component BA toward the flow axial which "pinches" and confines the charged aerosol; preventing the liquid droplets from reaching the duct wall; thus avoiding the wetting of the duct wall by the liquid, which might otherwise cause a short circuit.

The equation of the charged aerosol surface curve is derived similarly to the derivation for the duct curve (36) using a different value of j; but with the same orifice radius r_(b) at z=b; since the gas and the charged aerosol mix and flow through the same orifice; as shown in FIG. 2: Hence r_(za) =r_(b) (z/b)^(j/2) (40.1)

It is a feature of this invention to utilize a major portion of the duct curve (36) and the charged aerosol curve (40.1); and to discard the region of the duct curve (36) from (z/b)=(a/b) to (z/b)=0.20 and to discard the region of the charged aerosol curve from (z/b)=(a/b) to (z/b)=0.10. In these regions, the curves diverge rapidly to a large value. For a duct entrance diameter D_(a) =10 mm; an orifice diameter of 3.7 mm may be used. This results in an increased flow rate, and heat/electric power conversion of about 1 KW, as set forth hereinafter. The 10 mm dia. straight duct meets the duct curve (36) at about (z/b)=0.25, as shown in FIG. 1. A short curved surface joins the inlet duct surface and converging duct surface.

6.7 Uncoupled Flow

For uncoupled flow the continuity equations (30) apply only to the gas in the duct but not to the charged aerosol. By "uncoupled flow" it is meant that the area exponent for the duct J and the area exponent for the charged aerosol, j, in general, have different values; but the gas and the charged aerosol, assumed to have a charged droplet of near zero mobility, are both subject to the same velocity gradient. This velocity gradient has a maximum value of η≈1.5, as established by equation (11). With these conditions, a power ratio of 50 is obtained for χ=1 and 100 for χ=2, according to (24) and (25).

FIG. 2 shows plots of the curve of the duct surface, and the curve of the charged aerosol surface; using J=-1.41 for the duct and employing the continuity equation (30), and with j=-0.65 for the curve of the charged aerosol surface. The equations were derived for (b/a)=100:

    r.sub.d =1.85(z/b).sup.-0.705                              (41)

    r.sub.a =1.85(z/b).sup.-0.325                              (42)

The velocity ratio curve of FIG. (2) is plotted from

    (U.sub.z /U.sub.b)=(z/b).sup.1.5                           (43)

Referring to FIG. 2, power conversion starts at the maximum radius of the charged aerosol surface which is at (z/b)=0.08, near the entrance of the duct. The gas flow enters the converging duct at (z/b)=0.25, where it reaches the surface at a velocity ratio of 0.125. If the velocity at z=b is 86.5 m/s, then at (z/b)=0.25 the velocity.

6.8 Charged Aerosol Formation

Referring to FIGS. 1 and 2, a charged aerosol 10 is shown forming from a conductive liquid jet 7 and spreading within a cone to the entrance to the conversion space at z=a. An electric field is applied between the exciter electrode 8 and the liquid jet 7, by the exciter voltage source 9. As previously described⁶ the liquid jet enters the gas stream at a greater temperature than the gas, the electric field causes the jet to break into many small charged droplets charged by induction, the gas is locally heated, and the charged mutually repel each other; and a sudden expansion occurs, causing an increase of the velocity of the charged aerosol, momentarily resulting in a large value of the velocity gradient exponential, n, in the vicinity of the apex of the cone at the liquid jet. Consider the equation (14) for the E field ratio ζ at z=a (for example z=0.01 mm at the liquid jet for a 10 mm diameter exciter ring 8); (b/a)=100 now taken between the jet and the exciter electrode and j=2 for a conical expansion of the charged aerosol. Equation (14) for the E field ratio is independent of j at z=a. Evaluating (14), it is shown that for a current ratio of 1/ζ=50, n=1.5 is required; and this value of n may be obtained by the sudden expansion of the charged aerosol as it is formed, as above described. The maximum voltage generated in the cone χ=0.0125, calculated from (14) and (18), is very small, and is readily exceeded by the applied voltage between the exciter ring and the liquid jet 7, which results in large currents being drawn by the formation of the charged aerosol. Previous work has shown that the electric power expended in the formation of the charged aerosol is 0.1-1% of the output electric power.⁷

6.9 Isothermal Heat/Electric Power Conversion

The gas requires a time of about t=L/U_(ave) =10⁻² /50=2×10⁻⁴ sec to traverse the conversion space. However, with the small diameter of the charged droplets (about 0.3 μm) the time for temperature equilibrium, is much less, about 3.5×10⁻⁸ sec.⁸ Hence to assure that the droplets deliver their heat content to the gas gradually along the conversion space, at the rate required by electric power conversion, it is preferred to have the gas and the charged aerosol droplets enter the conversion space at the same temperature. In this manner, as electric power is converted from the heat content of the charged aerosol gas along the conversion space, this heat power is provided by (1) the enthalpy heat power decrease of the gas as it cools by a temperature ΔT; that is

    ΔQ.sub.1 =m.sub.g C.sub.p ΔT                   (44)

and by (2) the change i heat power content ΔQ₂ of liquid droplets:

    ΔQ.sub.2 =m.sub.o C.sub.o ΔT=αm.sub.g C.sub.o ΔT (45)

The temperature changes ΔT in (44) and (45) are equal and simultaneous. Hence the total heat power provided to the gas for the almost isothermal conversion to electric power is

    ΔQ=ΔQ.sub.1 +ΔQ.sub.2                    (46)

The gas and the charged liquid droplets enter the conversion space at z=a at temperature T_(a), and leave the conversion space at z=b at temperature T_(b).

Since ##EQU12## in which ○1 is the heat provided by the cooling liquid droplets, and ○2 is the heat provided by the cooling gas.

The heat ΔQ provided to the gas is converted to electric power and to the kinetic power of the gas as it increases velocity. The kinetic power component is small compared to the electric power output. The temperature of the gas and the charged droplets decrease, and the gas pressure and density decrease. This provides the electric power and the kinetic power outputs from the converter space.

Because the gas and liquid droplets are introduced at the entrance to the conversion space, at the same temperature, in temperature equilibrium, the nearly isothermal heat/electric power conversion is independent of the droplet diameter. The droplet diameter and electric charge per droplet have values such as to provide optimum power transduction, the charged droplet having a small mobility in the electric field.

The expansion is polytropic, not ideally isothermal. As previously described⁹ for a gas and a metal droplet which have a temperature decrease ΔT<<T, the expansion is approximately isothermal; from (46.2): ##EQU13##

Evaluating (47) for liquid metal tin droplet, and nitrogen gas, and assuming that the liquid droplet and gas heat contents are equal ##EQU14##

Putting these values into (47)

    1η(ρ.sub.1 /ρ.sub.2)=7(ΔT/T)             (49)

for

    (ρ.sub.1 /ρ.sub.2)=2.01;

    lη2.01=0.7

and

    (Δ

    T/T)=0.1                                                   (50)

Hence the heat/electric power conversion is approximately isothermal; thus for T=1800° K., ΔT=180° K. For smaller pressure ratios ΔT is smaller; thus for (ρ₁ /ρ₂)=1.5, (ΔT/T)=0.06 and ΔT=103° K.

It was previously shown¹⁰ that the electric power converted in a charged aerosol is:

    P.sub.e =1/2ε.sub.o b.sub.1.sup.2 (b.sub.g K.sub.a).sup.2 δ.sub.r.sup.2 UAΓ                             (51)

This equation may be expressed in terms of the mass flow rate, m_(g) with important consequences:

In (51) multiply and divide through by the density δ, and by R and make these substitutions: ##EQU15##

Define as a constant:

    K.sub.e =1/2ε.sub.o b.sub.1.sup.2 28.97/δ, 273×8314 (55)

    K.sub.e =1/2×8.854×10.sup.-12 ×(3×10.sup.6).sup.2 ×28.97/1.273×273×8314

    K.sub.e =4.00×10.sup.-4                              (56)

From (51)-(56) inclusive:

    P.sub.e =K.sub.e (b.sub.g K.sub.a).sup.2 m.sub.g RTδ.sub.r Γ=4×10.sup.-4 (b.sub.g K.sub.a).sup.2 m.sub.g RTδ.sub.r Γ                                                   (57)

The equation for electrothermodynamic power conversion is:

    P.sub.Q =m.sub.g RT ln (δ.sub.1 /δ.sub.2)=P.sub.e +P.sub.k (58)

In (58), P_(e) is the heat power converted to electric power; and P_(k) is the heat power converted to kinetic power, and exhausted from the orifice at (z/b)=1. The change in kinetic power is:

    P.sub.k =1/2(m.sub.o +m.sub.g)(U.sub.b.sup.2 -U.sub.a.sup.2) (59)

    P.sub.k =η.sub.k ·1/2m.sub.g (1+α)U.sub.b.sup.2 (59.1)

where η_(k) ≈1 since

    U.sub.b >>U.sub.a                                          (60)

The sonic velocity U_(s) is given by:

    U.sub.s.sup.2 =γRT                                   (61)

For a charged aerosol gas, the sonic velocity is U_(s) ¹

    1/2(1+α)m.sub.g U.sub.s.sup.12 =1/2m.sub.g U.sub.s.sup.2 (62)

Hence the sonic velocity of a charged aerosol gas is: ##EQU16## in terms of the Mach No M: ##EQU17##

From (59) to (64):

    P.sub.k =η.sub.k 1/2(1+α)m.sub.g ·M.sup.2 [1/(1+α)]γRT                                  (64.1)

    P.sub.k δ1/2m.sub.g RTγM.sup.2 ; since η.sub.k ≈1 (65)

From (57), (58) and (65):

    m.sub.g RT ln (δ.sub.1 /δ.sub.2)=K.sub.e (b.sub.g k.sub.a).sup.2 m.sub.g RTδ.sub.r Γ+1/2m.sub.g RTγM.sup.2 (66)

    ln (δ.sub.1 /δ.sub.2)=K.sub.e (b.sub.g k.sub.a).sup.2 δ.sub.r Γ+1/2γM.sup.2                   (66.1)

For P_(k) to be negligible:

    1/2γM.sup.2 << ln (δ.sub.1 /δ.sub.2)     (67)

if (δ₁ /δ₂)=2 and γ=1.4:

    1/2(1.4)M.sup.2 << ln 2=0.7

    M.sup.2 <<1                                                (67.1)

Example: ##EQU18## From (48) and (64) ##EQU19## The velocity at the entrance to the duct at (z/b)=0.2 is, from (43):

    U.sub.z =86.5(0.2).sup.1.5 =7.73 m/s                       (69)

Subject to condition (67) the last term in (66) may be neglected; whereupon this important criterion results:

    δ.sub.R Γ=[1/K.sub.e (b.sub.g k.sub.a).sup.2 ] ln (δ.sub.1 /δ.sub.2)                            (70)

Example: ##EQU20## The criterion (70) is satisfied for

    δ.sub.r =20; Γ=50; χ=1

or

    δ.sub.r =10, Γ=100 χ=2                     (72)

According to the criterion equation (70), the electric power output is dependent on the electric breakdown factor (b_(g) k_(a))² which preferably is 1 or more. For air at ambient temperature b_(g) k_(a) =1, but may be about 0.5 at 1800° K. As previously reported, certain additive gases will increase b_(g) k_(a), and may be employed in the practice of this invention. Alternatively for smaller b_(g) k_(a), a smaller pressure ratio will serve.

The increase in entropy due to electric power conversion is from (57): ##EQU21##

The mass flow rate m_(g) is controlled by the pressure ratio ρ₁ /ρ₂ and the orifice diameter D_(b) at z=b. The gas flow rate through the duct is subject to the continuity equation (41), which for a positive value of η≈1.5 in a converging flow with (ρ₁ /ρ₂)=1.5 within a duct results in a surface curve subject to the condition

    J=0.09-η,                                              (74)

and a radius r_(z) varying according to an exponent J/2 in (z/b)⁻⁰.705. This duct surface provides a velocity ratio at (z/b) of (U_(z) /U_(b))=(Z/b)^(1/)η required for maximum conversion of the heat power to electric power. The charged aerosol occupies a volume within the duct having a different surface; following a surface subject to the condition j+n=0.84; for which, from the equation (26), there is calculated a maximum power ratio of 50, for χ=1. Thus the duct surface establishes an accelerating flow, and the charged aerosol surface a charge distribution such as to result in a maximum power ratio.

Example

Given gas:air, electric breakdown factor

    b.sub.g k.sub.a =1 ambient

or

    0.5 at 1800° K.                                     (75)

Sonic velocity at 1800 K. ##EQU22##

Orifice velocity for α=4.76

    U.sub.b =0.2 U.sub.s '=0.2×850/(1+4.76).sup.1/2 =70.8 m/s (77)

Gas density at orifice ##EQU23## The mass flow rate takes into account the gas flow, which is decreased by the mass flow rate of the liquid tin. The area of orifice, for a diameter D_(a) =3.7×10⁻³ m=3.7 mm

    A.sub.b =(π/4)(3.7×10.sup.-3).sup.2 =1.075×10.sup.-5 m.sup.2 (79)

    m.sub.g =δ.sub.g A.sub.b U.sub.b =3.86×1.075×10.sup.-5 ×70.8=2.94×10.sup.-3 kg/s                     (80)

From (58) the electric power output is

    P.sub.e =P.sub.Q -P.sub.k =ηP.sub.Q =ηm.sub.g RT ln (ρ.sub.1 /ρ.sub.2)                                             (81)

for η=0.94 and (ρ₁ /ρ₂)=2.1; that is,

    P.sub.k =(1-η)P.sub.Q =0.06 P.sub.Q                    (81.1)

    P.sub.e =0.94×2.94×10.sup.-3 ×(8314/29)×1800 ln 2=1059 watts                                              (82)

Using the criterion (70) to find δ_(r) :

    δ.sub.r Γ.sup.1 =[1/K.sub.e (b.sub.g k.sub.a).sup.2 ] ln (ρ.sub.1 /ρ.sub.2)=[1/4×10.sup.-4 ×1] ln 2.1=1854 (83)

For Γ=50 (24)

    δ.sub.r =37                                          (84)

For Γ¹ =100 (25)

    δ.sub.r =18.5                                        (85)

These require different values of j: for χ=1 and χ=2, -0.635 and -0.85 respectively; which implies that the charged aerosol surface curves vary. Such variation may be obtained by increasing or decreasing the exciter voltage, varying the load resistor (16), or the pressure ratio ρ₁ /ρ₂.

The gas temperature increase decreases the electric breakdown factor b_(g) k_(a). Thus if b_(g) k_(a) =0.5 then the pressure ratio may be decreased; for example to 1.2 for the same value of δ_(r) ; but this will decrease the electric power output to 1/4×1059=264 watts. Hence to maximize electric power output and efficiency, there is an optimum operating temperature, and a trade-off, which is an engineering decision.

6.10 Liquid Volume Flow Rate

The liquid mass flow rate is by definition:

    m.sub.o =αm.sub.g                                    (86)

By continuity

    m.sub.o =δ.sub.o U.sub.o A.sub.o =δ.sub.o ν.sub.o =αm.sub.g                                           (87)

    ν.sub.o =(α/δ.sub.o)m.sub.g                 (88)

The mass flow rate of the gas is:

    m.sub.g =δ.sub.g ν.sub.g                          (89)

The gas density is:

    δ.sub.g =Ω.sub.r δ.sub.1 (273/T)δ.sub.r (90)

The ratio R_(og) of the liquid volume to the gas volume is:

    R.sub.og =ν.sub.o /ν.sub.g =(α/δ.sub.o)m.sub.g (m.sub.g /δ.sub.g)=α(δ.sub.g /δ.sub.o)     (91)

From (90) and (91):

    R.sub.og =[(δ.sub.1 /δ.sub.o)Ω.sub.r 273](αδ.sub.r /T)                             (92)

Example:

Find R_(og) for the following conditions

α=4.76 (from 48)

δ₁ =1.273 kg/m³ density of air at 1 atm and 273° K.

δ_(o) =5750 kg/m³ liquid tin

Ω_(r) =28/28.97 for N₂ relative to air rel. mol wt

δ_(r) =20

T=1800 K.

Solution: from (92)

    R.sub.og =[(1.273/5750)(28/28.97)273](4.76×20/1800)=3.08×10.sup.-3 (93)

That is, the volume ratio of liquid tin to N₂ is about 3/1000 but the mass ratio is 4.76.

6.11 Area and Diameter of the Liquid Jet Orifice

    A.sub.o =m.sub.o /δ.sub.o U.sub.o =αm.sub.g /δ.sub.o U (94)

Example: Given

    m.sub.g =2.94×10.sup.-3 Kg/S from                    (80)

    α=4.76 from                                          (48)

    U.sub.o =10 m/s from                                       (9)

    A.sub.o =4.76×2.94×10.sup.-3 /5750×10=2.4×10.sup.-7 m.sup.2                                                   (95)

6.12 Jet Power Input

    P.sub.o =1/2m.sub.o U.sub.o =1/2αm.sub.g U.sub.o.sup.2 =1/2×4.76×2.74×10.sup.-3 ×10.sup.2 (96)

    P.sub.o =0.70 watts                                        (97)

6.13 Droplet Diameter and Electric Charge

For optimum mobility¹¹. ##EQU24##

That is, there should be electron for each 100 Å Dia of the droplet.

    N=m.sub.o /(π/6)d.sub.o.sup.3 ·δ.sub.o =αm.sub.g /(π/6)d.sub.o.sup.3 δ.sub.o                      (99)

    I=N N.sub.e ·e=[αm.sub.g /(π/6)d.sub.o.sup.3 δ.sub.a ][10.sup.8 d.sub.o ]e                       (100)

    I=[(6/π)10.sup.8 e/δ.sub.a ][αm.sub.g /d.sub.o.sup.2 ]=5.32×10.sup.-15 (αm.sub.g /d.sub.o.sup.2)   (101)

    d.sub.o =7.29×10.sup.-8 (αm.sub.g /I).sup.1/2  (102)

Example: Given

    P.sub.e =1059 watts                                        (82)

    V=10.sup.5 volts

    I=P/V=1059/10.sup.5 =1.059×10.sup.-3 amps

    α=4.76                                               (48)

    m.sub.g =2.94×10.sup.-3 kg/s                         (80)

Find: The droplet dia. d_(o). From (102):

    d.sub.o =7.29×10.sup.-8 (4.76×2.94×10.sup.-3 /1.059×10.sup.-3).sup.1/2 =0.27×10.sup.-6 m

    d.sub.o =0.27 μm                                        (103)

From (98):

    N.sub.e =0.27×10.sup.-6 ×10.sup.8 =27 electrons/droplet (104)

7. The Generator

FIG. 3 shows a cross-sectional view of an electrothermodynamic generator of this invention full scale, comprising a 9 element array on the vertical axis. Each element generates about 1100 Watts; the total output being about 10 KW at 100,000 volts and 0.1 amp. D.C. The generator is one component in a Marks/Ericsson Cycle previously described.¹² It is maintained at its operating temperature T; for example 1800° K. within a heated enclosure, or furnace. The discussion relates to the generator, only.

The generator comprises two stainless steel end plates 21 and 22 which are at ground potential. An inlet pipe 23 supplies a gas; for example, air, nitrogen, or other suitable gas or mixture of gases¹³ through the entrance port 24 of the end plate 21 to a common chamber 25 in which the gas is temporarily held under a suitable supply pressure, for example 120 atmospheres and at T.=1800° K. The end plate 21 is welded to an intermediate plate 26, which contains a chamber 27 between the two plates. The chamber 27 receives the conductive liquid 28, which enters the chamber at a temperature T₁ ; for example 2400 K., and at a pump pressure of for example 130 atmospheres.

The end plates may be fastened together by screws to hold all the members including the insulating discs, 32, 33, and 34 under compression. This has not been illustrated herein for simplicity, but will be well understood by those skilled in engineering design.

The electrothermodynamic power conversion occurs in the ducts 30^(a), 30^(b), 30^(c), etc.; which contain charged aerosols 31^(a), 31^(b), 31^(c) etc located within the duct surfaces, separated therefrom, and operating as uncoupled flows, as particularly described in Section 6.7 herein. The ducts 30 are preferrably formed in a 10 mm thick sapphire plate 32. Sapphire is the best insulator, and non-reactive material suitable for these temperatures. For generators operating at lower temperatures other materials such as quartz may be used instead. Insulating spacer discs 33 and 34 may be joined to the plate 32; by suitable joints 37 and 38; and to the metal plates 22 and 26. Such joints have been described^(14'). The charged aerosols 31 are discharged by ions of opposite sign 4, emitted from points 39^(a), 39^(b), 39^(c), etc. which are mounted on an open mesh metal screen 40 held in place between elements 32 and 34. Since large voltages are involved the 100K lead 41 is passed through an insulating tube 42, which is sealed at 43 into a hole drilled in the side of plate 32. Similarly the exciter electrode 8, which is a conductive metal film shown in FIG. 2, on the entrance face of the insulating disc 32 is connected to the lead 44, which is brought out through the insulating tube 45, sealed to a second hole in disc 32, at 46.

After the charged aerosol is discharged by the points 39 the neutral droplets no longer repel each other, and coalesce into the liquid pool 47, which is withdrawn through the tube 48, as liquid flow 49, which is circulated through the liquid pump and heating coils back to the generator at 28, as previously described.¹⁵

The gas enters the generator at pressure p₁ and temperature T_(a) ≈T₁ through tube 23 in end plate 24; and exits through tube 50 as a gas stream, now separated from the liquid, at pressure p₂ and temperature T_(b). For example (p₁ /p₂)=1.5 or 2 and (T_(a) -T_(b))=ΔT 180 to 100K [ΔT/T₁ =(180/1800)=0.1] at 1800° K. Each duct may produce about 1100 watts (82); that is a total of about 10 WK for the 9 elements contained within the volume of a 7.5 cm cube.

In the practice of this invention, as described herein there may be employed a variety of gas and conductive liquids, each chosen for the selected ranges of temperature and pressure.

The gas may be selected from air, nitrogen, oxygen, hydrogen, helium and the other monoatomic gases; and compound gases such as the fluorocarbons, sulfur hexafluoride. Electron absorbing gases such as chlorine and carbon tetrachloride may be used to increase the electric breakdown strength of the gases, as previously described.¹³

The conductive liquids which may be employed include water, ethyl and the higher alcohols, various other organic fluids, and metals such as tin, mercury, gallium, metal compositions such as Woods Metal, Gallium-Niobium.

The objective is to obtain a liquid nonevaporating droplet in an inert gas having a higher breakdown strength, which is nonreactive with the materials, and stable at the maximum operating temperature chosen. With a high breakdown strength gas, the operating pressure is decreased.

The electrothermodynamic system with the equations and configurations (Section 6.9) hereinabove described, is substantially reversible, and may be advantageously employed by an ETD Compressor with charged water droplets in a Marks/Ericsson Cycle previously described.¹⁰

As above described, the gas and the charged aerosol are formed and enter the compressor at z=a, and discharge at z=b; with the velocity U increase from z=a to z=b with the velocity exponent n. However, in the compressor, the voltage is applied from z=a to z=b, such that V=V at z=a and V=0 at z=b, and the electric field intensity is reversed, so that electric power is supplied to the charged aerosol gas, to increase its pressure from p₂ to p₁, and increase its kinetic power. The gas-charged aerosol almost isothermally absorbs the electric power, and converts it to heat power to the charged water droplets, which serve as the heat-sink of the system at 300 K. The input electric power to the ETD Compressor supplies power to compress the charged aerosol gas from p₂ to p₁ while it is being cooled isothermally by the charged water droplets, and also to increase the velocity and kinetic power of the charged aerosol gas; which as in the ETD Generator is a small fraction of the electric power supplied to the compressor. After discharge, the liquid water is collected and sent to the water heat-sink heat exchanger.

The gas supplied to the ETD Compressor at pressure p₂ is supplied from the heat exchanger. Since its absolute temperature has decreased by the ratio 1800/300=6; the density is decreased by this factor. The isothermal work done on the gas is also decreased by this factor. Consequently if the ETD Generator has 12 modules producing 12 KW at 1800° K., the compressor then has only 2 modules, using 2 KW, and is very compact, in the same proportion.

This design is illustrative, and may be varied according to design specifications by those skilled in the art. Various modifications may be made herein without departing from the scope of this invention; for example, the D.C. output may be inverted to 60 HZ AC for general utility.

REFERENCES

1. U.S. Pat. No. 4,395,648 issued July 26, 1983 to Alvin M. Marks. For background see Cols. 2-8 inclusive and for the Bibliography see Cols 51 and 52

2. ibid: Method III Cols 31-44 inclusive

2.1 ETD Converter With Conical Jet CIP of 1. Copending Case Ser. No. 479,986 filed Mar. 29, 1983.

3. ibid: Col. 32 Equ (122) and Col. 34 Equ (145)

4. ibid: Col. 52 reference 7 also see Col. 4 lines 10-30

5. ibid: 4

6. ibid: Col. 44 lines 63-68

7. ibid: Bottom Col. 4, top Col. 5

8. ibid: Col. 48 Equ (220)

9. ibid: Cols. 46-47

10. ibid: Col. 34 Eq. (145)

11. ibid Col. 52 Refs. 7, 26

12. ibid Cols. 44, 45 FIGS. 26, 28-30 incl.; Cols. 46-51 incl.

13. ibid Col. 52 Ref. 7 Cols 21-22 Table No. 1.

14. ibid Col. 11 lines 38-48 inclusive Note: other materials may be plasma sprayed, and a higher temperature braze metal may be employed for the 1800-2400 K. range

15. ibid Col. 46, FIG. 28. 

Having fully described my invention, what I wish to claim is:
 1. In an electrothermodynamic generator for the substantially isothermal conversion of the internal heat power of a flowing gas charged aerosol to electric power, a source of gas, said gas entering said generator at a temperature T₁ and pressure ρ₁ a gas source, a duct connected to said gas source, said duct having a flow axis in the Z direction, said duct converging from an entrance diameter D_(a) at Z=a to an exit orifice in said duct having a diameter D_(b) at Z=b, a source of conductive liquid at a pressure and temperature exceeding that of the said gas, a minor orifice for said liquid located on said axis upstream from said duct entrance, said conductive liquid issuing as a jet from said minor orifice, an exciter electrode located at the said duct entrance, an exciter electric potential source, said potential source being connected across said conductive liquid and said exciter electrode to cause a charged aerosol droplets to form from said liquid jet between said minor orifice and said exciter electrode, a discharge electrode located downstream from the said exit orifice, which injects ions of opposite sign into said charged aerosol, the heat content of said charged aerosol droplets being at least equal to the heat content of the gas in which it is contained, constituting an internal heat source to maintain the charged aerosol gas at nearly constant temperature, said converging duct causing the velocity and kinetic power of said gas-charged aerosol to increase as it moves along the flow axis from Z=a to Z=b; said gas decreasing in pressure from ρ₁ at Z=a to ρ₂ at Z=b, the ratio of temperature decrease ΔT to the absolute temperature T₁ being less than about 0.1, whereby the charged aerosol gas almost isothermally converts its heat power to electric power, said kinetic power being a small fraction of said electric power; a load resistor, said load resistor being connected between said discharge electrode and said conductive liquid source, said electric power being delivered to said load.
 2. An electrothermodynamic generator according to claim 1, in which said charged aerosol lies on an inner surface inside the gas flowing in the duct, and separated therefrom, except in the vicinity of the exit orifice, a velocity gradient established by said gas flowing in said duct establishing the said velocity gradient being the same for the said gas and said charged aerosol, but the spatial distribution of the electric charge density of said aerosol being determined by its said inner surface, whereby a current and power outputs are increased at a given operating voltage.
 3. An electrothermodynamic generator according to claim 2 in which the velocity exponent n has a value of about 1.5, and in which the area exponent of the duct is about J=-1.41 governed by continuity, and in which the area exponent of the charged aerosol is about j=-0.65.
 4. An electrothermodynamic generator according to claim 2 in which the curve of the duct surface has the equation:

    r.sub.d =1.85(z/b).sup.-0.705

and in which the curve of the charged aerosol surface has the equation:

    r.sub.a =1.85(z/b).sup.-0.325.


5. An electrothermodynamic generator according to claim 1 in which the said gas flowing in said duct is electrically uncoupled said charged aerosol flow being separated from the walls of said duct.
 6. An electrothermodynamic generator according to claim 1, an insulating disc and an array of said converging ducts formed in said insulating disc, whereby a multiplicity of said generators are connected in parallel to increase the electric power output.
 7. An electrothermodynamic generator according to claim 1, said generator being incorporated in a Marks/Ericsson Cycle.
 8. In an electrothermodynamic generator according to claim 1, a conductive liquid chosen from the class consisting of water, alcohol and a liquid metal including tin, Woods Metal, mercury, gallium, gallium-noibium.
 9. In an electrothermodynamic generator according to claim 1, a gas selected from the class consisting of air, oxygen, nitrogen, hydrogen, and other monoatomic gases.
 10. In an electrothermodynamic generator according to claim 9, said gas containing additives to increase electric breakdown strength selected from the class consisting of fluorocarbons, sulfur hexafluoride, chlorine, carbon tetrachloride.
 11. In an electrothermodynamic generator according to claim 2, in which the said charged aerosol has a spatial charge distribution along said flow axis providing an electric potential along the said flow axis having a substantially decreased electric field intensity everywhere along said axis, thereby enabling an increase in the electric current and electric power output by several orders of magnitude, with substantially no increase in the output voltage.
 12. A method for the substantially isothermal electrothermodynamic conversion of the internal heat-power of a flowing charged aerosol gas to electric power, which comprises the steps of providing a gas at a temperature T_(a) and pressure ρ_(a), confining said flowing charged aerosol gas within a converging duct, said duct having a flow axis ZZ', providing a conductive liquid at a pressure temperature exceeding that of said gas, causing said liquid to flow through a minor orifice to form a conductive liquid jet upstream from said duct entrance, an exciter electrode at the entrance to said duct, applying an electric field between said exciter electrode and said jet, forming charged aerosol droplets in said flowing gas stream, said charged aerosol droplets having a mass flow rate M_(g), said flowing gas having a mass flow ratio α, the ratio α=m_(o) /m_(g) being such that the heat content of the said liquid droplets at least exceeds that of the said gas; the heat contents of said liquid droplets and the said gas constituting an internal heat source, passing said charged aerosol through the said converging duct from Z=a to Z=b on said axis respectively at the entrance and the orifice at the exit of the said duct, causing the said gas to decrease in temperature and pressure from T_(a) and ρ_(a) at Z=a to T_(b) and ρ_(b) at Z=b respectively, the value of ΔT/T being less than about 0.1 to obtain substantially isothermal operation, where ΔT=T_(a) -T_(b) ; increasing the velocity and kinetic power of the gas, and converting the difference between said heat power and kinetic power to electrical power, discharging said charged aerosol at the orifice of said duct by spraying said charged aerosol with ions of opposite polarity from a collector electrode, said kinetic power being a small proportion of said electric power, delivering said electric power to an electrical load.
 13. A method for the substantially isothermal electrothermodynamic conversion of the internal heat power of a flowing charged aerosol gas according to claim 12, in which the charged aerosol flows along a surface within the duct surface which have a common exit diameter at z=b; in which the flowing gas establishes a common increasing velocity space gradient, and in which the charged aerosol surface establishes the charge distribution to obtain a large power ratio of the order of 50 to
 100. 14. As an article of manufacture, an electrothermodynamic power converter comprising a first end plate, an inner chamber in said end plate, a gas inlet connection, and liquid inlet connection in said first end plate, a plurality of minor orifices in said end plate communicating with said inner chamber and said liquid inlet, a central insulating plate containing a plurality of converging ducts in an array, first and second insulating spacers, an exciter electrode comprising a conductive metal film as the first surface of the said insulating plate, a plurality of discharge points on a supporting screen mounted between said second surface of said insulating space and said second spacer, a second end plate, liquid and gas outlet connection in said second plate, said first end plate, first insulating spacer, insulating plate, second insulating spacer, and second end plate being mounted and sealed together to form a pressure-containing vessel; a first major chamber between said first end plate and said first surface of said insulating plate to receive the gas and provide space for the formation of the charged aerosol from the plurality of liquid jets from said minor orifices, and a second major chamber between the second surface of said insulating plate and said second end plate, to provide space to receive the discharged aerodol droplet from the gas, and to provide a liquid pool to receive and coalesce the said discharged liquid droplets, and to exhaust the gas through its outlet; a first insulated tube containing the high voltage outlet lead, sealed to the side of said plate, said outlet lead being connected to said discharge electrodes; and a second insulated tube containing the high voltage exciter lead, said exciter lead being connected to the conducting layer on the inner face of said insulating plate.
 15. An article of manufacture according to claim 14 in which said end plates compress the said insulating plate and spacers.
 16. An article of manufacture according to claim 14, said article being incorporated in a Marks/Ericsson Cycle.
 17. In a Marks/Ericsson Cycle according to claim 7,an electrothermodynamic compressor for the substantially isothermal compression of a flowing charged aerosol-gas at a temperature T₂ from pressure ρ₂ to pressure ρ₁, a gas source, a duct connected to said gas source, a source of electric power, said duct having a flow axis in the z direction, said duct converging from an entrance diameter D_(a) at z=a, to an exit diameter in said duct having a diameter D_(b) at z=b, a source of conductive liquid at a pressure exceeding that of the gas, a minor orifice for said liquid located on said axis upstream from said duct entrance, said conductive liquid issuing as a jet from said minor orifice, an exciter electrode located at the said duct entrance, an exciter electric potential source, said potential source being connected across said conductive liquid and said exciter electrode to cause charged aerosol droplets to form from said liquid jet between said minor orifice and said exciter electrode, the heat content of said charged aerosol droplets being at least equal to the heat content of the gas in which it is contained, constituting an internal heat-sink to maintain the charged aerosol gas at a nearly constant temperature, said converging duct causing the velocity and kinetic power of said charged aerosol gas to increase as it moves along the flow axis from z=a to z=b, said gas increasing in pressure from ρ₂ at z=a to ρ₁ at z=b, the ratio of temperature increase ΔT to the absolute temperature T₂ being less than about 0.1, a discharge electrode at said duct exit to discharge said aerosol by ions of opposite electric sign thereto, said electric power source being connected between said conductive liquid and said discharge electrode, whereby said charged aerosol almost isothermally converts the said electric power to said increase in pressure and kinetic power of said flowing gas, a heat-sink heat exchanger, said discharged droplets being collected and supplied to said heat-sink heat exchanger.
 18. In an ETD compressor according to claim 17, charged aerosol droplets comprising water liquid. 